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ABSTRACT 

We present temporal and spectral characteristics of X-ray flares observed from 
six late-type G-K active dwarfs (V368 Cep, XI Boo, IM Vir, V471 Tau, CC 
Eri and EP Eri) using data from observations with the XMM-Newton obser- 
vatory. All the stars were found to be flaring frequently and altogether a total 
of seventeen flares were detected above the "quiescent" state X-ray emission 
which varied from 0.5 to 8.3 xlO^^ ergs s~^. The largest flare was observed in 
a low activity dwarf XI Boo with a decay time of 10 ks and ratio of peak flare 
luminosity to "quiescent" state luminosity of 2. We have studied the spec- 
tral changes during the flares by using colour-colour diagram and by detailed 
spectral analysis during the temporal evolution of the flares. The exponential 
decay of the X-ray light curves, and time evolution of the plasma temperature 
and emission measure are similar to those observed in compact solar flares. 
We have derived the semiloop lengths of flares based on the hydrodynamic 
flare model. The size of the flaring loops is found to be less than the stel- 
lar radius. The hydrodynamic flare decay analysis indicates the presence of 
sustained heating during the decay of most flares. 

Key words: X-ray:stars - stars:activity - starsxoronae - stars:flare - starsdate- 
type 



1 INTRODUCTION 

Flares are events in which a large amount of energy is released in a short interval of time. 
Such events take place at almost all frequencies of the electromagnetic spectrum and have 
been observed in the Sun as well as in many types of cool stars (Garcia- Alvarez, et al. 
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2002). These stellar flares could be radiating several orders of magnitude more energy than 
a solar flare. Optical flares in UV Get (dMe) type stars are a common phenomenon. The 
flares produced by other stellar sources (e.g. RS CVn and BY Dra) are usually detected 
only in the UV or X-rays. These UV and X-ray flares show extreme luminosities and very 
hot temperatures (> lOMK). Even though, flares in these stars present many analogies with 
the solar flares, there are also signiflcant differences, such as the amount of energy released. 
A general model for stellar flares has emerged from numerous solar flare studies (Giidel 
2004). The flare reconnection region, located somewhere at large coronal heights, primarily 
accelerates electrons and ions upto MeV energy (Dennis & Schwartz 1989). The accelerated 
electrons participate along the magnetic flelds into chromosphere where they heat the cool 
plasma to coronal flare temperatures, thus evaporating a part of the chromosphere into 
the corona. Early X-ray observations with the EXOSAT observatory of flare stars revealed 
examples of two different type of flares: 1) Impulsive flares which are like compact solar 
flares, and 2) Long decay flares which are like two-ribbon solar flares (Pallavicini, Tagliaferri 
& Stella 1990). The compact flares are less energetic (~ 10^° ergs s~^), short in duration 
(< 1 h), and conflned to a single loop while the long-decay flares are more energetic (~ 
10^^ ergs s~^), of long duration (^ 1 h), and release the energy in an entire arcade of loops 
(Garcia- Alvarez 2000). The most hkely flare process relates to an opening up of magnetic 
flelds and subsequent relaxation by closing the open field lines. 

Analysis of light curves during fiares can provide us with insights into the characteristics 
of the coronal structures and, therefore, of the magnetic field (e.g., Schmitt & Favata 1999; 
Favata, Micela & Reale 2000a; Reale et al. 2004). Even though stellar flares are spatially 
unresolved, a great deal of information on the coronal heating and on the plasma structure 
morphology can be inferred from a detailed modeling of stellar flares; for instance, if sufficient 
data are available for moderately time-resolved spectral analysis, a study of the complete 
evolution of a flare can allow us to (a) infer whether the flare occurs in closed coronal 
structures (loops), (b) determine the size of the flaring structures, (c) determine whether 
continuous heating is present throughout the flare, and (d) put constraints on the location 
and distribution of the heating (see Reale et al. 2004). 

Active solar-type (G-K) dwarfs are known to possess magnetic flelds on their surfaces, 
and the fleld strengths are as large as several kilogauss, i.e. stronger than the strongest 
fleld observed on the Sun (e.g., Saar 1990; Johns-Krull & Valenti 1996). Gyclic behavior has 
been identifled in several late- type active stars (Baliunas et al. 1995, Olah & Strassmeier 
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2002), suggesting that dynamos similar to the solar dynamos are also operative in solar- 
type stars. These solar- type stars are less active than the dMe stars (so called flare stars). 
However, they show more frequent flaring activity than the late-type evolved stars (e.g. RS 
CVns; Giidel 2004). In this paper, we present analysis of archival data obtained from XMM- 
Newton observations of six G-K dwarfs, namely V268 Cep, XI Boo, IM Vir, V471 Tau, 
EP Eri and CC Eri with the aim to understand the spectral and temporal characteristics 
of X-ray flares in them. In subsequent papers we will examine the characteristics of X-ray 
flares in the subgiants and the giants. 

The paper is organized as follows: In ^ we give the basic parameters of stars in our 
sample, in ^we describe the observational data sets that we have analysed and the methods 
of data reduction, §1] contains our analysis and results, and in [|5] we present our discussion 
and conclusions. 

2 BASIC PARAMETERS OF SAMPLE STARS 

The basic parameters of stars in our sample are given in Table 1. V368 Cep (= HD 220140; 
V = 7.54 mag) is a G9V spectral-type chromospherically active star with photometric period 
of 2.74 d (Chugainov, Lovkaya & Petrov 1991, Chugainov, Petrov & Lovkaya 1993). XI Boo 
(= HD 131156; V = 4.55 mag) is a nearby visual binary, comprising a primary G8 dwarf 
and a secondary K4 dwarf with an orbital period of 151 yrs (Hoffleit & Jaschek 1982). In 
terms of the outer atmospheric emission, UV and X-ray observations show that the primary 
dominates entirely over the secondary (Hartman et al. 1979, Ayres, Marstad & Linsky 1981, 
Schimtt 1997, Laming & Drake 1999). XI Boo A (primary) is a slow rotator with rotation 
period of 6.2 d (Gray et al. 1996). A cool main sequence eclipsing binary (G5V+KV), IM Vir 
(= HD 111487; V = 9.69 mag ) has an orbital period of 1.31d (Malkov et al. 2006). V471 Tau 
(= BD-l-16 516; V = 9.71) is a rapidly rotating detached eclipsing binary first observed by 
Nelson & Young (1970). The system is composed of a cool main-sequence chromospherically 
active K2 V star (Guinan & Sion 1984) and a degenerate hot white dwarf separated by 3.4 
stellar radii. CC Eri (= HD 16157, V = 8.8) is a spectroscopic binary with orbital period 
of 1.561 d. It consists of a K7Ve primary and a dM4 secondary, with mass ratio of ~ 2 
(Strassmeier et al. 1993). The primary star corotates with orbital motion due to the tidal 
lock, and is one of the fastest rotating K dwarfs in the solar neighborhood. EP Eri (=HD 
17925, V = 6.0 mag) is very nearby (10.3 pc), very young (high Li I abundances), active K2- 
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type dwarf with a rotation period of 6.85 days (Cayrel de Strobel & Cayrel 1989; Cutispoto 
1992; Henry, Fekcl & Hall 1995). The presence of an unresolved companion in this star has 
been suggested by Henry et al. (1995) based on the variable widths of the photospheric 
absorption hues reported in the literature {vsini ranges from 3 to 8 km s~^; see Fekel 1997). 

3 OBSERVATIONS AND DATA REDUCTION 

The late-type dwarf stars in our sample were observed with the XMM-Newton satellite 
using varying detector set-ups. The XMM-Newton sateUite is composed of three co-aligned 
X-ray telescopes (Jansen et al. 2001) which observe a source simultaneously, accumulating 
photons in three CCD-based instruments: the twin MOS 1 and MOS 2 and the PN (Turner 
et al. 2001; Striider et al. 2001), all three detectors constituting the EPIC (European Photon 
Imaging Camera) camera. The EPIC instrument consists of three CCD cameras with two 
different types of CCD, two MOS and one PN, providing the imaging and spectroscopy in 
the energy range from 0.15 to 15 keV with a good angular (PSF = 6 arcsec (FWHM)) and 
a moderate spectral resolution {E/AE 20 — 50). Exposure time for each star was in the 
range of 30-60 ks. A log of observations is provided in Table 2. 

The data were reduced with standard XMM-Newton Science Analysis System (SAS) 
software, version 7.0 with updated cahbration files (Ehle et al. 2004). The preliminary pro- 
cessing of raw EPIC Observation Data Files was done using the epchain and emchain tasks 
which allow calibration both in energy and astrometry of the events registered in each CCD 
chip and to combine them in a single data file for MOS and PN detectors. The background 
contribution is particularly relevant at high energies where coronal sources have very little 
flux and are often undetectable. Therefore, for further analysis we have selected the energy 
range between 0.3 to 10.0 keV. Event list file was extracted using the SAS task evselect. The 
epatplot task was used for checking the existence of pile-up affecting the inner region of all 
stars. Only the star CC Eri was found to be affected by pile-up. X-ray light curves and spec- 
tra of all target stars were generated from on-source counts obtained from circular regions 
with a radius ~ 40-50" around each source. However, for the star CC Eri, X-ray hght curves 
and spectra were extracted with events taken from an annulus of 32"with inner radius of 
12"to avoid pile-up effect in the inner region. The background was taken from several source 
free regions on the detectors at nearly the same offset as the source and surrounding the 
source. 
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4 ANALYSIS AND RESULTS 

4.1 Multi-band X-ray light curves and characterization of flares 

The background subtracted X-ray light curves of six stars viz. V356 Cep, IM Vir, V471 Tau, 
CC Eri and EP Eri, as observed with the MOS and PN detectors are shown in Figures [H^a) 
to[T](f), respectively. All light curves are in the energy band 0.3-10.0 keV. The MOS and PN 
light curves are represented by solid and open circles, respectively. The temporal binning of 
light curves is 200 s for all sources. The PN light curve of the star IM Vir was affected by 
proton flares twice during the observations, flrst for 800 s after 10 ks from the beginning of 
the observation, and second at the end, after 40 ks of the observations. However, the MOS 
light curve was affected only at the end of the observation. For the star EP Eri, the high 
proton flare background region in the light curve is for a duration of 3 ks after 14.7 ks from 
the beginning of the observation. We have removed data during the proton flare background 
from the MOS and PN light curves of both the stars IM Vir and EP Eri. The star XI Boo 
was not observed by the PN detector, and only the light curve obtained with the MOS is 
presented here. 

The light curves of all the individual stars show variability on a time scale of ks, most of 
which resembles flaring activity. The customary deflnition of a flare is a significant increase in 
intensity, after which the initial or quiescent level of intensity is reached again. Such patterns 
are observed in the light curves of the sample stars shown in Figures [T](a) to[T](f), and where 
the "flare regions" are represented by arrows and marked by Fi, where i=l,2...17 refers to 
the flare number. The mean quiescent level count rates of the sources are taken from regions 
that are free from flares and are marked by Q in the Figure [H The start time, end time, 
flare duration and flux during quiescent and at flare peak for all the 17 flares observed are 
listed in Table 3. The peak and quiescent state count rates as given in Table 3 are converted 
into flux by using the WebPIMM^, where we assumed the plasma temperature of 1 keV. 
Flares F6, F7, F9, F14 and F16 were found to be long lasting 2.3 h) flares. Infact, the 
flare F7 of the XI Boo is longest flare observed with total duration of ~ 3.1 h. The peak flux 
in these long duration flares was found to be more than two times than the quiescent state. 
Other detected flares were shorter than 1.3 h. The shortest duration flare was detected one 
in V368 Cep (F3) and other in CC Eri (F12). In these short duration flares, peak X-ray flux 
was 1.3 - 1.8 times more than that of the quiescent state. 



^ http:/ /heasarc. gsfc.nasa.gov/Tools/w3pimms. html 
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To characterize the flares, we have fitted the hght curves of the flares with an exponential 
function 

c{t) = Aoe~'~^ + q (1) 

where c(t) is the count rate as a function of time t, to is the time of peak count rate, q 
is the count rate in the quiescent state, is the decay time of the flare, and Aq is the count 
rate at flare peak. The best fitted parameters for the flares are given in Table 4 for both the 
MOS and the PN data. Five flares, F6 and F7 of XI Boo, F9 of IM Vir, F14 of CC Eri and 
F16 of EP Eri the e-fold decay time was found to be more than 1 h. Flare F7 of XI Boo was 
found to be the longest decay flare (t^ = 10 ks) observed among all the seventeen flares. The 
flare F16 shows a very different structure when compared to the other flares, i.e. a sharp 
rise and a highly structured decay as if more flares were present during the decay. Though 
it was tough to locate the peak of flare F16 because of its complex structure, the decay time 
was determined to be 4 ks. For the remaining eleven flares the decay time was found to be 
in the range of 8 - 47 minutes. The decay phase of the flare F17 was not observed. The rise 
times of all flares were found to be less than 1 h. The rise and decay time of the shortest 
flare F3 of V368 Cep was found to be 0.5 ks. In XI Boo, after the flare F6 an active level U 
was identified, where average flux was 1.8 times more than that of the Q state (see Figure 
[D^b)), that was higher than the peak of the flare F5 indicating a very high activity level. 
Two similar highly active regions Ul and U2 were identified between the flares F13 and 
F14. The average fluxes during the regions Ul and U2 were 1.3 and 1.4 times more than 
that of the quiescent state. Presence of strong substructures in these active regions indicates 
superposition of a large number of flaring regions. 

To investigate the behavior of the flares detected in different energy bands, the light 
curves of the stars obtained with the MOS and PN are divided into three energy bands 
namely soft (0.3 - 0.8 keV), medium (0.8 - 1.6 keV), and hard (1.6 - 10 keV). The boundaries 
of the selected energy bands are chosen as the line free regions of the low resolution PN 
spectra. The hardness ratio HRl and IIR2 are defined by the ratio of medium to soft band, 
and hard to medium band count rate, respectively. The soft, medium, hard band intensity 
curves, and the hardness ratio curves HRl and IIR2 as a function of time are shown in the 
sub-panels running from top to bottom in left(MOS) and right(PN) of Figure [2] (a) for the 
star V368 Cep. Similar plots of intensity and hardness ratios observed from the stars XI Boo, 
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IM Vir, V471 Tau, CC Eri and EP Eri are shown in the Figures [2]^b) to[2]^f), respectively. 
The hght curves in individual band also show significant variability in all sources. We have 
determined the e-folding decay time of observed fiares in the three bands by fitting equation 
([1]). The decay times, rise times, count rates at fiare peaks, and quiescent state count rates 
at soft, medium and hard band for PN data are given in the Table 5. For most of the fiares, 
e-folding decay time in the soft band was found to be more than the decay time in the 
medium and the hard band. However, for the fiare F3 of V368 Cep and F7 of XI Boo, it 
appears that decay time in the medium band was more than in the soft band. For the fiares 
F6, F7 and F9 the decay time in the soft band was found to be more than 3.0 ks than that 
of the hard band. However, for other fiare events the decay time in the soft band was < 2 ks 
than that in the hard band. Two smaller fiare like structures were seen just after the peak 
of fiare F13 of CC Eri, in the hard band but were not observed in the soft and medium band 
light curves (see Figure [2]^e)). The rise time for all fiares in all band was found to be well 
within la level (Table 5). 

The fiare peak to the quiescent state fiux ratio in the hard band was found in between 
2-16, which is more than that of the similar ratios in the medium and the soft band (see 
also Table 5) for all fiares. In terms of the peak fiux, the fiares F6 and F7 of XI Boo, and 
F13 and F14 of CC Eri were stronger in the medium band than in the soft and the hard 
band. The variation in the hardness ratio are indicative of changes in coronal temperature. 
An X-ray fiare is well defined by a rise in the temperature and subsequent decay. For the 
fiares F6, F7, F9, F12, F13 and F14, both HRl and HR2 varied in the similar fashion to 
their light curves (see Figure [2]). This implies that an increase in the temperature at fiare 
peak and subsequent cooling. The structures observed during the decay phase of the fiares 
F2 of V368 Cep and F13 of CC Eri in the hard band were also observed in the HR2 curve. 
This could probably be due to emergence of fiares during the decay phase of these fiares, 
which were comparatively stronger in the hard band. During the fiares F2, F3 and F4 of 
V368 Cep, and F16 of EP Eri the HRl did not show any variation but HR2 varied according 
to the fiare intensity. However, HRl and HR2 did not show any significant variation during 
the fiares Fl, F5, F8, FIO, Fll, F15 and F17. The PN data have more count rates than the 
MOS and have similar characteristics to that of MOS, therefore, for further analysis in this 
paper we use only the PN data. 



8 J. C. Pandey and K. P. Singh 
4.2 Colour-colour diagrams 



Plots of hardness ratio in the form of colour-colour (CC) diagram can reveal information 
about spectral variations and serve as a guide for a more detailed spectral analysis. For 
this purpose, the light curves and the hardness ratio curves of the flares were divided into 
different time segments covering their rising and decaying phases as shown in Figure [1], 
where Ri(i=l,2..) represents the rising phase, Di(i=l,2,..) represents the decay phase and 
Q represents the quiescent state. We determined the HRl and HR2 from the PN data for 
each of these segments and plotted the values in Figure [31 To understand the observed 
behavior of HRl and HR2 in terms of simple spectral models, we generated the soft(0.3- 
0.8 keV), medium(0.8-1.6 keV) and hard(l. 6-10.0 keV) band count rate using the "fakeit" 
provision in the XSPEC, using the most recent Canned Response Matrices downloaded 
from http://xmm.vilspa.esa.es/external/xmm_SW_cal/calib/epic_files.html, and simple two 
temperature plasma models. The two temperature APEC models were mixed such that the 
emission measure were in the ratio EM2/EM1 of 0.4, 0.6, 0.8, 1.0 and 1.2 (based on spectral 
results in §4.3 and §4.4). The abundances (0.25) and hydrogen column density (10^° cm~^) 
were kept fixed for each set of the EM2/EM1 (see §4.31 and §4.4p . Further by keeping kTi 
fixed at different values ranging from 0.2 to 1.2 keV in steps of 0.2 and varying kT2 from 
0.4 to 2.4 keV such that kT2 > kTi, we generated count rates in the soft, medium and hard 
bands. These predicted count rates were used to generate the hardness ratios. The families 
of curve thus generated are over-plotted on the data points in Figure [3] to understand the 
CC diagram. The kT2 in Figure [3] increases from bottom to top of each curve. 



Figure [3] (a) shows the plot between HRl and HR2 for a set of EM2/EM1 = 1.0 for the 
star V368 Cep. In the case of the flare Fl, 'Dl' intersects the generated CC curves for which 
kTi is 0.6 keV, however, 'D4' intersects the curve generated for kTi = 0.2 keV. A similar 
pattern is seen for the other flares F2 and F4 of the stars V368 Cep, where the top of the 
flare is located near the CC curve for a high temperature component whereas the colours 
near the end of the flare intersect the CC curve corresponding to a lower temperature. The 
flare F2 shows a very high temperature during its rise phase 'Rl' and the decay phase 'D5', 
as both are located between the CC curves of 1.0 and 1.2 keV. The decay phase D6 of the 
flare F2 is located near the 0.8 keV CC curve. The peak seen during the decay phase D6 
of flare F2 (Fig[3](a)) can be explained by a rise in the temperature of coronal plasma. It 
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appears that the quiescent state 'Q' of the star V368 Cep has a cool temperature between 
0.2 to 0.4 keV. 

The peak of the largest flare 'F9' observed in the star IM Vir appears to have a tem- 
perature between 1.0 - 1.2 keV. As the flare decayed the temperature was also observed to 
decrease. A similar trend is seen during the flares observed in the stars XI Boo (Figure [3](b)), 
V471 Tau (Figure EKd)), CC Eri (Figure [S^e)), and EP Eri (Figure EKf)). The observed CC 
curves for the star IM Vir, V471 Tau and CC Eri were well matched with the generated fam- 
ily of CC curves for a set of EMi/ EM2 = 1.0. However, the observed CC diagram of EP Eri 
were well matched with generated CC curves for a set of EM1/EM2 = 0.4. It appears that 
high temperatures are needed during the flare peak, and the plasma temperatures decrease 
as flare decays. To, further confirm and quantify the spectral changes during the flare, we 
have also performed a detailed spectral analysis of the flares and the quiescent states of the 
stars. The results of these analyses are presented below. 

4.3 Quiescent state X-ray spectra 

X-ray spectra for each star during the quiescent state were analysed using XSPEC version 
12.3 (Arnaud 1996). Data were cleaned of proton flares by removing the affected time periods 
for the stars IM Vir and EP Eri. Spectral analysis of EPIC data was performed in the energy 
band between 0.3-10.0 keV for the star V368 Cep, XI Boo, V471 Tau and CC Eri, while for 
the stars IM Vir and EP Eri spectral analysis was performed between 0.3-5.1 keV. Individual 
spectra were binned so as to have a minimum of 20 counts per energy bin. The EPIC spectra 
of the stars were fitted with a single (IT) and two (2T) temperature coUisional plasma model 
known as APEC (Astrophysical Plasma Emission Code, Smith et al. 2001), with variable 
elemental abundances. Abundances for all the elements in the APEC were varied together. 
The interstellar hydrogen column density (Nh) was left free to vary. For all the stars, no 
IT or 2T plasma models with solar photospheric (Anders & Grevesse 1989) abundances 
could fit the data, as unacceptably large values of xl were obtained. Acceptable APEC 2T 
fits were achieved only when the abundances were allowed to depart from the solar values. 
The best-fit 2T plasma models with sub-solar abundances along with the significance of 
the residuals in terms of Ax^ are shown in Figures Hl^a) to (f) for all the stars. Table 6 
summarizes the best-fit values obtained for the various parameters along with the minimum 
x'i = x"^/^ (where u is degrees of freedom) and the 90% confidence error bars estimated from 
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the minimum + 2.71. The derived values Nh from the spectral analysis were found to be 
lower than that of the total galactic HI column density (Dicky & Lockman 1990) towards 
the direction of that star. The cool and hot temperatures at the quiescent state of these 
stars were found in the range of 0.2 - 0.5 and 0.6 - 1.0 keV, respectively. The corresponding 
EM2/EM1 during the quiescent state of each star was found in the range of 0.8 - 1.2, except 
for the star EP Eri where it was found to be ~ 0.4. 

4.4 Spectral evolution of X-ray flares 

In order to trace the spectral changes during the flares, we have analysed the spectra of 
the different time intervals shown in Figure [TJ To study the flare emission only, we have 
performed 1-T spectral fits of the data, with the quiescent emission taken into account by 
including its best-fit 2-T model as a frozen background contribution. This is equivalent to 
consider the flare emission subtracted of the quiescent level, allows us to derive one "effective" 
temperature and one emission measure of the flaring plasma. The abundances were kept fix 
to that of the quiescent emission. The best fit spectral parameters for each flare segments 
are listed in Table 7. Figures [5] - [10] show the temporal evolution of the temperature and 
corresponding emission measure of the flares detected in the targeted stars, where '0' time 
corresponds to the flare peak time. 

4.4.1 V368 Cep 

Spectral data of flares Fl, F2 and F4 were binned into four segments. However, for the 
flare F3 data could be collected over only one segment. Both the temperature and the 
corresponding emission measure were found to decrease from the flare peak to the quiescent 
state (see also Table 7). As shown in Figure [5](a), the highest temperature for the flare 
Fl was 14 MK and decreased along the decay path to 6.8 MK. A similar trend was also 
found in the emission measure, where it increased by a factor of seven (see Figure [5] (b)). 
The temperature and emission measure for the flare F2 were peaked simultaneously during 
the decay phase D5. However, for the flare F4 both parameters were peaked during the rise 
phase R2. The maximum luminosity was found during the flares F2 and peaked at 1.18 x 10'^° 
ergs s~^, which is ~ 1.7 times more than that of the quiescent state. 
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4.4.2 XI Boo 

To see the flare evolution, flare F6 and F7 have binned into five and six intervals, respectively. 
For the flare segments R2, D5, D6, D7 and D8, the IT APEC fit gave high value of xl (=1-5 
for R2, 1.75 for D5, 1.61 for D6, 1.55 for D7 and 1.34 for D8). The spectral fit to the data 
for these time segments improves the by using the 2T plasma model (see Table 7). The 
temperature and the corresponding emission measure of the cool component for these time 
segments were found to be constant at 0.6 ±0.1 keV and 4.0 ± 1.0 x 10^^ cm~^, respectively. 
The spectral parameters for hot component of the best fit 2T plasma model are reported in 
Table 7. The evolution of the temperature of hot component and the corresponding emission 
measure are shown in Figures [6] (a) and [6] (b), respectively. For the flares F6 and F7, both 
the temperature and corresponding emission measure were reached to the maximum value 
during the decay phase. As mentioned in Table 6 and 7, during the flares F6 and F7, XI 
Boo was ~ 2 times more X-ray luminous than that of the quiescent state. 

4.4.3 IM Vir 

The spectral data of the flares F8 and F9 were binned in two and four different time segments, 
respectively (see Figure[T]^c)) and spectra of each segment were fitted with IT APEC model. 
The plots of temperature and the corresponding emission measure are shown in Figures!?] (a), 
and (b), respectively. For the flare F9, the maximum fit temperature, was reached during the 
rising phase R2, somewhat earlier than the maximum of its emission measure. The emission 
measure was increased by a factor of four during the flare F9. The luminosity at the peak 
of flare F9 was found to be ~ 2 times more than that of the quiescent state. 

4.4.4 V47I Tau 

Two flares detected in the star V471 Tau were examined by analysing four different time 
intervals during the flares. Figures. [8] (a) and (b) show the plot of temperature and the 
corresponding emission measure during the flares FIO and Fll. For the flare FIO the tem- 
perature reached its maximum value of 17.8 MK and decreased along the decay phase (see 
Table 7). The corresponding emission measure was also peaked during the decay phase of 
the flare FIO (see Figure [8]^b) and Table 7). For the flare Fll, both the temperature and 
the emission measure follow similar trend to the flare light curve. The X-ray luminosity was 
increased by a factor of ~ 1.6 for the flare FIO and ~ 1.4 for the flare Fll. 
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4.4.5 CC En 

The evolution of temperature and the corresponding emission measure for three flares F12, 
F13 and F14 are shown in Figures [9] (a) and (b), respectively. Small variations were found 
in the temperature and the emission measure during the flares F13 and F14. Both the 
temperatures and the corresponding emission measures were found to decrease during the 
decay phase for all flares (see also Table 7). The temperature and the corresponding emission 
measure for all the flares were peaked during the decay phase. However, it appears that 
temperature is peaked before the emission measure for the flare F13. The peak luminosity 
during flares F13 and F14 was found to be 6.0 x 10^*^ ergs s~^ more than that of the quiescent 
state. 

4.4.6 EP En 

Spectra of different time segments of the flares are fitted with IT plasma models. Figures 
[To] (a) and (b) show the evolution of temperature and the corresponding emission measure 
for the flares F15 and F16. For the flare F16, both the temperature and the corresponding 
emission measure were peaked during the rise phase R2, remained constant upto the decay 
phase D2 and then decreased along the flare decay without any appreciable changes in the 
luminosity. The luminosity during both flares was found to be similar to that of the quiescent 
state (see Table 6 and Table 7). 

4.5 Loop modeling of the X-ray flares 

Flares can not be resolved spatially on a star. However, by an analogy with solar flares and 
using flare loop models, it is possible to infer the physical size and morphology of the loop 
structures involved in a stellar flare. The most widely used but approximate methods for 
analyzing stellar flares are: i) Two-ribbon flare method (Kopp & Poletto 1984), ii) Quasi- 
Static Cooling method (van den Oord & Mewe 1989), iii) Pure radiation cooling method 
(Pallavicini et al. 1990), iv) Rise and decay time method (Hawley et al. 1995), and v) 
Hydrodynamic method (Reale et al. 1997). The Two- ribbon flare model assumes that the 
flare decay is entirely driven by heating released by magnetic reconnection of higher and 
higher loops and neglects completely the effect of plasma cooling. The other four methods 
are instead based on the cooling of plasma confined in a single fiaring loop. In the hypothesis 
of fiares occurring inside closed coronal structures, the decay time of the X-ray emission 
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roughly scales as the plasma cooling time. In turn, the cooling time scales with the length 
of the structure which confines the plasma: the longer the decay, the larger is the structure 
(e.g. Haisch et al.l983). A loop thermodynamic decay time has been derived (van den Oord 
& Mewe 1989; Serio et al. 1991) as: 

V J? 

where Lg and T7 arc the loop half-length and the maximum temperature (T,nax) of the 
flaring plasma, in units of 10^ cm and 10^ K, respectively. The timescale above is derived 
under the hypothesis of impulsive heat released at the beginning of a flare. 

The hydrodynamic model includes both plasma cooling and the effect of heating during 
flare decay. Reale et al. (1997) presented a method to infer the geometrical size and other 
relevant physical parameters of the flaring loops, based on the decay time and on evolution of 
temperature and the emission measure (EM) during the flare decay. The empirical formula 
for the estimation of the unresolved flaring loop length has been derived as: 



^'^^ /(C) > 1 (3) 



120/(0 

where Td is the decay time of the light curve, and f{C) is a non-dimensional correction 
factor larger than one. ^ is the slope of the decay path in the density-temperature diagram 
(Sylwester et al. 1993) and is maximum (~ 2) if heating is negligible and minimum (~ 0.5) 
if heating dominates the decay. The correction factor f{Q and loop maximum temperature 
are calibrated for the XMM- Newton EPIC spectral response and are given as (Reale 2007): 



^1 

/(C) = + 1-36 (for 0.35 < C ^ 1-6) (4) 



Trna. = 0.13Tj,f (5) 

Other physical properties of the flaring plasma can be inferred from the analysis of flare 
data. The analysis of X-ray spectra provides values of the temperature and emission measure. 
From the emission measure (EM) and the plasma density (rig), the volume of the flaring loop 
is estimated as (Reale 2002): 
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V-'-^ (6) 



where V is the loop volume. It has been shown that in the equilibrium condition the loop 
scaling laws hold, linking maximum pressure, temperature, loop length and heating rate per 
unit volume as (Rosner, Tucker & Vaiana 1978, Kuin & Martens 1982): 

T^a. = 1.4 X l{}\pLf/^ (7) 



Eh ~ IQ-^T^l^L-"^ ergs s'^cm-'^ (8) 

The minimum magnetic field necessary to confine the fiaring plasma can be simply estimated 

as 



B = ^/Mp (9) 

Because of limited statistics, we have modeled only ten fiares. The results of the model 
parameters are listed in Table 8. In this table, column 1 is name of the star and the cor- 
responding fiare, column 2 is slope of density-temperature diagram ((^), column 3 contains 
loop maximum temperature based on spectral fit and equation (5), semiloop length is de- 
termined using the equation (3) and given in column 4, pressure in the loop at the fiare 
peak is estimated using the loop scaling law (see equation (7); column 5), column 6 contains 
maximum electron density at the fiare loop (we have assumed a totally ionized hydrogen 
plasma i.e. p = 2nekTmax', since pressure derived from equation (7) is a maximum value, 
therefore, a upper limit on density is estimated), volume of the fiaring plasma is estimated 
from maximum emission measure and from the estimated maximum electron density, and 
therefore, is treated as lower limit (equation (6); column 7), column 8 is heating rate per unit 
volume and estimated using equation (8), the minimum magnetic filed necessary to confine 
the fiaring plasma at loop apex B (see equation [9]) is given in column 9, the loop aspect 
ratio (3 {=r/L, r is radius of the loop cross section, assuming circular) estimated from the 
loop scaling law V = 2TTf3'^L^ and number of loops filling the fiare volume estimated assum- 
ing that (3 = 0.1 for a single loop are given in column 9 and 10, respectively. The model 
parameters for the different fiares from each star are presented and discussed below. 
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4.5.1 Flares Fl, F2 and F4 of V368 Cep 

The average temperatures of the loops, usually lower than the real loop maximum tempera- 
tures, are found from the spectral analysis of the data as given in §4.4[ According to equation 
([5]), the loop maximum temperatures for the flare Fl, F2 and F4 are found to be 25.5 ± 2.5, 
62 ± 8 and 40.3 ± 1.8 MK, respectively. Figure [11] (a) shows the density-temperature (n-T) 
diagram, where EM^/^ has been used as a proxy of density. In Figure [11] (a) the flares Fl 
,F2 and F4 are represented by solid circles , solid triangles and solid squares, respectively. 
The solid lines represent the best linear fit to the corresponding data, providing the slope 
These values of ( indicate the presence of a sustained heating during the decay of the 
flares Fl and F2. However, sustained heating is negligible during the decay of the flare F4. 
As discussed in §4.5[ the e-folding time of the light curve and the slope of n-T path can be 
used to estimate the loop half-length. According to equation (3) the semiloop lengths for the 
flares Fl, F2 and F4 are determined to be 6.0±0.6x 10^ (= O.IR^), 19.8 ±8.6 (= O.AR^) and 
11.7 ± 1.3 X 10^ (= 0.2i?^) cm, respectively. These loop lengths are much smaller than the 
pressure scale height^ hp > 1.2 x 10^^ cm. Maximum pressure in the loop at the flare peak 
is estimated to be > 10'^ dyne cm~^ for the flares Fl, F2 and F4, respectively (see equation 
(7)). Assuming that the hydrogen plasma is totally ionized {p = 2nekTmax), the maximum 
plasma density in the loop at the flare peak is estimated in the order of 10^^ cm~^ for the 
flares Fl, F2 and F4, respectively. We computed a volume of ~ 1.78 x 10^° cm~^ for the flare 
Fl using the observed peak EM of 3.7 x 10^^ cm~'^. Similarly, loop volume is determined to 
be 7.2 X 10^^ cm^ for flare F2 and 7.4 x 10^^ cm^ for the flare F4. A hint for the heat pulse 
intensity comes from the flare maximum temperature. By applying the loop scaling laws 
and loop maximum temperature (see equation ([8|)) the heat pulse intensity for the flares 
Fl, F2 and F4 would be the order of 2.36, 4.75 and 3.06 ergs cm~'^ s~^, respectively. From 
the pressure of the flare plasma the magnetic field required to confine the plasma should be 
more than the 160 Gauss for all the flares. 

4.5.2 Flares F6 and Fl of XI Boo 

The evolution of flare F6 and F7 in n - T plane is shown in Figure [TT] (b), together with a 
least square fit to the decay phase. The resulting best-fit slopes {() for the decaying phase of 

^ defined as hp = 2kT/^g ~ 5000 X Tmax /{g/SQ), where T is plasma temperature in the loop, is the molecular weight of 
the plasmas and g is the surface gravity of the star. 
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flare F6 and F7 indicate that the flare F6 is driven by the time scale of the heating process, 
whereas sustained heating is neghgible during the decay of the flare F7. The intrinsic flare 
peak temperatures for flare F6 and F7 are, applying equation (5) to the observed maximum 
temperature , 18.9±0.8 and 25.3±2.7 MK, respectively. For the flare F6 the semiloop length 
is estimated to be less than the pressure scale height, hp = 10^^ cm. Maximum pressure and 
density in the loop at the peak of the flare F6 are estimated to be 130 ± 6 dyne cm~^ and 
2.5±1.0 X 10^'^ cm^^, respectively. Using equation[6l volume of the flaring plasma is estimated 
to be 1.1 X 10^^ cm'^. A minimum magnetic field of ~ 58 Gauss is required to confine the 
flaring plasma at the peak of the flare. Although, the value of ( for the flare F7 is outside 
the domain of the validity of the method (see Real 2007), the estimated loop parameters 
obtained from the hydrodynamic loop modeling are listed in Table 8. The estimated semiloop 
length for the flare F7 is smaller than the pressure scale height (=1.3 x 10^^ cm). 

4.5.3 Flare F9 of IM Vir 

Using equation ([5]) the maximum loop temperature for the flare F9 is found to be 46.6 ±3.7 
MK. Figure [TT] (c) shows the n-T diagram for the flare F9. From the best linear fit, we 
found the slope of the n-T diagram implies that the heating during the flare decay was as 
strong as in the flares Fl, F2 and F6. The resulting semiloop length is found to be ten 
times smaller than the pressure scale height {hp = 2.3 x 10^^ cm), which implies a loop in 
hydrostatic equilibrium with a plasma pressure of 2020 ± 1100 dyne cm~^ and maximum 
heating rate of ~ 2.07 ergs cm^'^ s^^. The maximum plasma density and loop volume for 
the flare are estimated to be 1.6 ± 0.8 x 10^^ cm~^ and 4.55 x 10^° cm^, respectively. The 
minimum confining magnetic field is estimated to be ~ 226 Gauss. 

4.5.4 Flare FIO of V471 Tau 

The maximum temperature in the loop at the flare peak was found to be 33.5 ± 3.3 MK for 
the flare FIO. (see Figure [TlTd) for the n-T diagram of this flare). The slope of the decay 
path of n-T diagram indicates the presence of sustained heating. Using equation (3) the 
semiloop length of the flare FIO is calculated to be ~ 0.4i?^ and is ~ 5 times smaller than 
the pressure scale height hp = 1.2 x 10^^ cm. Using these parameters and the scaling laws for 
the loops, the maximum pressure, the maximum density and the minimum magnetic field 
confining the flaring plasma in the loop at the flare peak, the volume of flaring plasma and 
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the heating rate per unit volume for FIO are estimated to be p = 630 ± 200 dyne cm~^, 
He = 6.8 ± 1.0 X 10^° cm~^ and B = 126 Gauss, 8.9 x 10^° cm^ and 0.46 ergs s~^ cm~^. 

4.5.5 Flares F13 and FI4 of CC Eri 

Using equation ([5]) the maximum temperature in the loop responsible for the flares F13 and 
F14 are found to be 25.7 ± 2.9 and 27.0 ± 2.1 MK, respectively. Figure [TT] (e) shows the 
n-T diagram for these flares. The slopes in n-T diagram for the flares F13 and F14 indicate 
the presence of strong sustained heating, and showing that the observed decay is driven by 
the time-evolution of the heating process. The value of C at the low extreme of the error 
bar for the flare F13 is compatible with the lower asymptotic value for which equation (4) 
can be applied. Therefore, an upper limit for semi-length of the loop can be derived by 
using the value of C, at the high extreme of the error bar. The resulting semi-length for the 
flares F13 and F14 are < 1.4 x 10^°(= 0.3/0.5i?J and 2.0 ± 1.0 x W^{= 0.4/0. 7/?J cm, 
respectively. These values of the loop length are much smaller than the pressure scale height 
hp = 1.1 X 10^^ for flare F13 and hp = 7.3 x 10^° cm for flare F14. A magnetic field more 
than 95 Gauss is required to confine the flaring plasma at the top of both flares. A detailed 
analysis of these flares is also given in Crespo-Chacon et al. (2007). The estimated loop 
parameters for both flares are found well within la level to that of Crespo-Chacon et al. 
(2007). 



4.5.6 Flare Fl 6 of EP Eri 

The trend in n-T diagram suggests the presence of significant heat during the decay of 
the flare F16 (See Figure [11] (f) and Table 8). The temperature at the peak of the flare 
is determined to be 20.2 ± 2.2 MK. The derived semiloop length L = 2.4 ± 0.5 x 10^° 
cm (= 0.4i?j,) is less than the pressure scale height hp = 9 x 10^^ cm. Using the same 
approach as discussed above, the maximum pressure and the maximum density are derived 
at p = 130 ± 50 dyne cm^^ and Ue = 2.3 ± 0.5 x 10^° cm^^. Using the maximum emission 
measure, the loop volume is estimated to be 5.38 x 10'^*' cm'^. The confining magnetic filed 
is estimated as S > 57 Gauss. 
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5 DISCUSSION AND CONCLUSIONS 



We have carried out an analysis of XMM-Newton observations of six G-K dwarfs. Light 
curves revealed flaring in all the stars on time scale ranging from seconds to kiloseconds, 
with various peak strengths. A total of 17 flares were identified in these G-K dwarfs. The 
decay time of these flares ranges from 0.5 to 10 ks. A similar range of flare decay times in the 
pre-main-sequence stars of Pleiades cluster was found by Stclzcr, Neuhauser & Hambaryan 
(2000). Among the 17 flares presented here, only four flares (F6 and F7 of XI Boo, F9 of IM 
Vir and F14 of CC Eri) appear to be long decay flare(T(i ^ 1 hr). Infact, the flare F7 of XI 
Boo was one of the longest duration flares observed with ~ 10 ks. Such long decay flares 
have so far been reported in CF Tuc (r^ = 79.2 ks; Kriister & Schmitt 1996), EV Lac (r^ = 
37.8 ks; Schmitt 1994), Algol (r^ = 30.2 ks; Ottmann & Schmitt 1996), AD Leo (r^ = 7.9 
ks; Favata, Micela & Reale 2000a) and some pre-main sequence stars (Stclzcr, Neuhauser & 
Hambaryan 2000, Favata et al. 2005). This flare classification is purely based on the decay 
time. It has been shown that even much longer flares can occur in a single loop (Favata et al. 
2005), while relatively short flares with heating dominated decay are probably arcade flares. 
For the flares F6, F9 and F14, the decay path is driven by the sustained heating, therefore, 
these flares can not be classified as long decay flares. However, the sustained heating during 
the decay of the flare F7 is negligible. Thus, these flares can be classiflcd as arcade. The flare 
decay time for the remaining 12 flares are similar to that of the solar compact flares. The 
morphological differences in these different types of flares indicate the different processes 
of energy released. In compact flares energy is probably released only during an impulsive 
phase, whereas in the two-ribbon flares a prolonged energy release is apparently required to 
explain their long decay time (Pallavicini, Serio & Vaiana 1977; Pallavicini et al. 1988; Priest 
1981; Poletto, Pallavicini & Kopp 1988). The rise time of these flares has been found to be 
less than 1 ks, which is similar to the rise time of impulsive flares observed in the M dwarfs 
(Pallavicini et al. 1990). However, for some flares (F6, F7 and F9) the rise time was found to 
be more than 1 ks. We found that the decay times of the flares observed here are more in the 
soft band than in the medium band or the hard band. Some flares, e.g., F2, actually showed 
a hard peak during the decay. But even then the in soft band being more than in the 
hard band. This could probably be due to the softening of the spectrum during the decay 
due to the plasma cooling i.e. emission gradually exists form high energy band and enters 
more deeply in the soft energy band or the higher energy electrons streaming down from 



A study of X-ray flares - I. Active late-type dwarfs 19 

the coronal heights to the chromosphere that can heat the plasma to higher temperature, 
lose their energy faster than the lower energy electrons. Most of the flares observed in our 
sample were found to have more peak flux in the soft band than in the medium and hard 
band. 

The total energy released during the X-ray flares observed in the sample of G-K dwarfs 
is observed to be in the range of 2.3 x 10^^ to 6.1 x 10^^ ergs. This shows that these flares are 
ten times more energetic than the flares in the M-type dwarfs (Pallavicini et al. 1990), and 
10^ to 10^ times more energetic than the solar flares (Moore et al. 1980, Wu et al. 1986). 
However, they are not as energetic as a large flare observed in Algol (Efot = 10^^ ergs; White 
et al. 1986), but are as energetic as the flares observed in G9 dwarf ZS 76 (Efot = 1-4 x 10^^ 
ergs; Pillitteri et al. 2005). The large flare, F7, of XI Boo was the largest in terms of both 
the decay time and the ratio of the peak luminosity to the quiescent luminosity (a factor of 
2.2). However, given the low quiescent emission (Lx = 5.1 x 10^^ ergs s~^) from the XI Boo, 
this flare was not particularly prominent in terms of either peak X-ray luminosity (Lx = 
1.1 X 10^^ ergs s~^) or the total X-ray energy {Efot = 6.4 x 10^^ ergs). Other large flare F9 of 
IM Vir shows a large peak X-ray luminosity (1.6 x 10'^° ergs s~^) as well as the total energy 
released {Etot = 6.12 x 10^^ ergs). 

Inspite of the fact that the peak flare luminosity may vary by a factor of 2 (see, e.g., IM 
Vir in Table 7) , there is a good correlation between the peak flare luminosities (Lx/) and the 
quiescent state stellar X-ray luminosities (Lxg). The form of the correlation is log Lx/ = 1.04 
log Lxg + 1.3, with the correlation coefficient, r, of 0.9 (see Figure fT2l) . A similar correlation 
was found by Pallavicini et al. (1990) for the M dwarfs. This indicates that there is a direct 
connection between the flaring and the quiescent X-ray emission, which is probably due to 
their common emission mechanism for X-rays. 

Flares F2-F3-F4, F6-F7 and FlO-Fll show the similar structure i.e. before ending the first 
flare the next flare starts. Similar loop systems have been observed to the flare on the Sun 
(e.g. so-called Bastille-day flare; Aschwanden & Alexander 2001), and in a stellar analogue 
dMe star Proxima Centauri (Reale et al. 2004). In these two events, a double ignition in 
nearby loops was observed or suggested, and the delay between the ignitions appears to scale 
with the loop sizes. Similar type of double exponential decay was also found in the young 
stars ZS 76 (Pillitteri et al. 2005) and Cygnus 0B2 (Albacete Colombo et al. 2007). This 
implies that in these dwarfs a steady active region undergoes a strong magnetic reconnection 
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event resulting in an intense flare, and is followed by an arcade of reconnected loops that 
slowly decay. 

We have performed time resolved spectroscopy of these flares. The coronal spectrum 
during the flare can be represented with a IT model, with quiescent state taken into account 
as a frozen background contribution. During the flares Fl, F9 and Fll the emission measure 
was increased by a factor of five or more. Such large variation in emission measures were also 
seen during the flares detected in AB Dor (Maggio et al. 2000). Both the temperature and 
emission measure show the well-defined trends i.e. the changes in the temperature and the 
emission measure are correlated with the variations observed in the light curves during the 
flares. Similar trends were also seen during the flares detected in the pre-main sequence stars 
of the Orion nebula cluster (Favata et al. 2005) and star-forming complex L1551 in Taurus 
(Giardino et al. 2006). Reale, Peres & Orlando (2001) found that the height of emission 
measure distributions are variable during the different phases of the solar flares, while its 
width and the peak temperature of the distribution undergo much smaller changes. The 
peak flare temperature is found in the range of 20 - 60 MK for all the flares. These values 
are intermediate between those of flares observed on active stars like Algol {Tmax ~ 100-150 
MK, Favata et al. 2000b) and AB Dor {T^ax ~ 140-170 MK, Maggio et al. 2000) and those 
found in flares on the dMe star AD Leo (range 20-50 MK, Reale & Micela 1998; Favata et al. 
2000a), and on solar- type Pleiades (20-40 MK, Briggs & Pye 2003), and young stellar objects 
{Tmax ~ 80-270 MK; Favata et al. 2001,2005). For most of flares both emission measure and 
temperature peaked simultaneously. However, for few flares (F9, Fll, F13), it appears that 
the temperature was evolved before the emission measure. Similar delay is often observed 
both in solar flare (Sylwester et al. 1993) and in flares from the stars Algol(van Den Oord 
& Mewe 1989), EV Lac (Favata et al. 2000a), AB Dor (Maggio et al. 2000), and YY Gem 
(Stelzer et al. 2002). This is probably due to an impulsive flare event , in which loop does 
not reach equilibrium conditions, the density begins to decay later than the temperature. 

We have modeled the flares using the hydrodynamic model (see §4.5j Reale et al. 1997) 
based on the decay phase of the flare. The derived semiloop lengths for all flares are found 
to be in the range of 0.6 — 8.0 x 10^° cm. Alternatively, Reale (2007) derive the semiloop 
length from the rise phase and peak phase of the flare as: 

L', ^ (10) 

where Lg is semiloop length in the unit of 10^ cm, ip = Tmax/Tu, T^ax is maximum tem- 
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perature in the unit of 10^ K, Tm is temperature at density maximum and tM is time in 
the unit of 10^ s at which density maximum occurs. The time of the maximum emission 
measure is a good proxy for Tm- In case where both temperature and emission measure were 
simultaneously peaked during the rise phase, only rise time was used in the estimation of 
the semiloop length. The semiloop lengths estimated from this approach for the flares F4 
{L'q = 13.6 ± 3.2 cm), ¥6^ = 16.8 ± 1.8 cm), F13(L(, = 6.6 ± 1.4 cm), FU{L'q = 11.4 ± 3.1 
cm) and F16(L9 = 20.8 ± 5.8 cm) are found to be consistent with that of estimated from 
the decay phase analysis (see also Table 8). The estimated loop length, L'g, for the flares F2 
(9 ±2 cm), F7(36±8 cm) and FIO (7 ±2) are found to be less than that of determined from 
decay phase. However, for the flare F9 of IM Vir the estimated loop length (Lg = 41 ± 13) 
is two times more than that of estimated from decay phase. The inconsistency in the esti- 
mation of the loop length from the two approaches is probably due to the involvement of 
other different coronal loops during the decay of the flare or the heat pulse triggering the 
flare is not a top-hat function. This new analysis allows us to derive the loop lengths of the 
flares for which either time resolved spectral information is not available or decay phase is 
not observed. These flares are Fll of V471 Tau, and F15 and F17 of EP Eri. The estimated 
semiloop lengths for these flares Fll, F15 and F17 are 6.2 x 10^°, 2.1 x 10^° cm and 3.8 x 10^ 
cm, respectively. Therefore, the magnetic structures confining the plasma for all the observed 
flares in G-K dwarfs are smaller than the star themselves, and are not as large as observed in 
a single giant HR 9024 (L = i?^/2; Testa et al. 2007) and some pre- main- sequence analogues 
(L >> R^; Favata et al. 2005, Giardino et al. 2006). 

For all the flares, the estimated maximum electron density under assumption of a totally 
ionised hydrogen plasma is found in the order of 10^°"^^ cm~'^. This is compatible with the 
values expected for the plasma in coronal condition (Landini et al. 1986). To satisfy the 
energy balance relation for the flaring as a whole, the maximum X-ray luminosity must be 
lower than the total energy rate {H = Eh-V; see §4.5p at the flare peak. The rest of the input 
energy is used for thermal conduction, kinetic energy and radiation at lower frequencies. For 
the flares F2 and F4 of V368 Cep the maximum X-ray luminosity observed is about 34% 
of H. Similarly, for the flares F7 of XI Boo and FIO of V471 Tau the maximum X-ray 
luminosity is about 30% of H. However, for other flares (see Table 7) only 16 - 22 % of if 
is observed as a peak luminosity. These values are in agreement with those reported for the 
solar flares, where the soft X-ray radiation only accounts for upto 20% of total energy (Wu 
et al. 1986). In comparison, the fraction of X-ray radiation to the total energy has been 
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found to be 15% and 35% for the M dwarfs Proxima Centauri and EV Lac, respectively 
(Reale et al. 2004; Favata et al. 2000c). Applying loop scaling law V = 27r/?^L^, and if the 
detected flares are produced by a single loop, their aspect ratio {(3) were estimated in the 
range of 0.1 to 0.3 for the flares F2, F4, F7, F9, F13, F14 and F16 (see Table 8). Similar 
cross section was also observed for the solar coronal loops, for which typical values of (3 are 
in the range of 0.1 -0.3. If we assume (3 = 0.1, these flares occurred in the arcades, which are 
composed of 2 to 11 loops (see also Table 8). However, for the flares F6 and FIO the value 
of P (assuming single loop) is estimated to be 0.51 and 0.37, respectively. In comparison to 
the solar coronal loops, the cross section for these flare loops are very large. It appears that 
the flares F6 and FIO occurred in the arcades that contain 26 and 14 loops, respectively. 

Therefore, we conclude that the observed flares in G-K dwarfs are similar to the solar 
arcade flares, which are as strong as M dwarfs and are much smaller than the flare observed 
in dMe star, giants and pre-main-sequence analogous. 
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Table 1. General pro]:)erties of stars in the sample 



Name 


Spectral Type 


V 


Period 


Distance 






(mag) 


(days) 


(pc) 


V368 CEP 


G9V 


7.54 


2.74 


19.7 


XI Boo 


G8V/K4V 


4.55 


6.2 


6.6 


IM VIR 


G5V/KV 


9.00 


1.3086 


60.0 


V471 Tau 


K2V/dW 


9.71 


0.52 


46.7 


CC Eri 


K7Ve/M4V 


8.76 


1.56 


11.5 


EP ERI 


K2V 


6.00 


6.85 


10.3 
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Tabic 2. Log of oljservations with the XMM-Newton 



Star 
Name 


Instrument 
(mode,filter) 


Start time 
date (UT) 


Exposure time 
(s) 


offset (') 


V368 CEP 


MOSl(SW,thick) 


2003-12-27(19:23:43) 


29762 


0.045 




M0S2(SW,thick) 


2003-12-27(19:23:43) 


29768 






PN(SW, thick) 


2003-12-27(19:28:55) 


29570 




XI Boo* 


M0S1(T, medium) 


2001-01-19(11:25:06) 


58397 


0.150 




M0S2(SW, medium) 


2001-01-19(11:25:06) 


58997 




IM VIR 


M0S1(FF, medium) 


2004-07-15(04:49:45) 


53488 


9.093 




M0S2(FF, medium) 


2004-07-15(04:49:45) 


53487 






PN(FF, medium) 


2004-07-15(05:32:47) 


51077 




V471 Tau 


M0S1(LW, medium) 


2004-08-01(06:52:10) 


60678 


1.052 




M0S2(LW, medium) 


2004-08-01(06:52:10) 


60682 






PN(SW, medium) 


2004-08-01(06:57:35) 


60471 




CC Eri 


M0S1(PFW, thick) 


2003-08-08(08:40:36) 


39463 


0.065 




M0S2(T, thick) 


2003-08-08(08:40:27) 


39212 






PN(PFW, thick) 


2003-08-08(09:21:36) 


36703 




EP Eri 


M0S1(SW, thick) 


2004-01-20(17:53:18) 


50013 


0.015 




M0S2(SW, thick) 


2004-01-20(17:53:19) 


50013 






PN(SW, thick) 


2004-01-20(17:58:29) 


50013 





*XI Boo was observed for 895 s by EPIC PN detector 

Note: FF: Full frame, FPW: FuUPrime window PPW: PrimePartial Window, SW: 
Small window, T: timing, LW: Large window 
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Table 3. Time interval of 17 flares observed in six G-K dwarfs. 



Object Flare Time-interval Duration Peak flare flux^ Quiescent state flux^ 

(ks after the observation start) (ks) 10~^^ ergs s~^ cm~^ lO"^'^ ergs s~^ cm~^ 



V368 Cep 


Fl 


.. - 4.0 


4.0 


2.28 


1.30 




F2 


19.2 - 22.5 


3.3 


2.39 






F3 


22.5 - 23.9 


1.4 


2.36 






F4 


23.9 - 27.8 


3.9 


1.90 




XI Boo 


hr, 


8.1 - 12.7 


4.6 


2.87 


2.14 




F6 


12.7 - 21.1 


8.4 


5.66 






F7 


21.1 - 32.4 


11.3 


5.57 




IM Vir 


F8 


4.5 - 9.3 


4.8 


0.20 


0.13 




F9 


9.3 - 18.9 


9.6 


0.25 




V471 Tau 


FIO 


14.4 - 18.5 


4.1 


0.39 


0.24 




Fll 


18.5 - 23.0 


4.5 


0.31 




CC Eri 


F12 


5.8 - 7.5 


1.7 


2.38 


1.83 




F13 


8.3 - 12.5 


4.2 


3.45 






F14 


27.1 - 36.5 


9.4 


3.43 




EP Eri 


F15 


4.5 - 8.5 


4.0 


0.62 


0.48 




F16 


35.1 - 44.2 


8.9 


0.65 






F17 


47.4 - .. 




0.66 





count rates are converted into the flux using WebPIMMS (http://heasarc.gsfc.nasa.gov/Tools/w3pimms.html) 



Table 4. Parameters obtained from fitting of equation 1 to the flare light curves obtained 



with MOS and PN detectors in the energy band of 0.3-10.0 keV 



Object 


FN 






MOS 








PN 




Lx" 


Etot 


Name 




Ao 






Tr 


q 


Ao 






Tr 


q 


(10^3 


(10^^ ergs) 






(counts s 






(s) 


(counts s~^) 


(counts s 




(3) 


(s) 


(counts s~^) 


ergs s~^) 




V368 Cep 


Fl 


7.0 ± 0, 


.2 


1213 ± 74 




3.6 


9.0 ± 0. 


.2 


1226 ± 86 




5.1 


31.3 


12.5 




F2 


6.8 ±0, 


.1 


1658 ± 141 


316 ± 43 




9.1 ± 0, 


.3 


2387 ± 341 


347 ± 63 




39.6 


13.1 




F3 


6.5 ± 0, 


.2 


451 ± 49 


479 ± 44 




8.7 ±0. 


.2 


496 ± 36 


616 ± 79 




9.6 


9.6 




F4 


4.5 ± 0, 


.3 


1117 ± 110 


649 ± 181 




7.0 ± 0. 


.2 


1232 ± 133 


681 ± 91 




29.9 


11.7 


XI Boo 


FB 


3.42 ± 0, 


.12 


2829 ± 950 


879 ± 213 


2.66 












0.6 


0.23 




F6 


6.9 ± 0, 


.1 


4393 ± 333 


3169 ± 269 














4.3 


3.61 




F7 


6.78 ±0, 


.06 


9885 ± 362 


3361 ± 244 














6.7 


6.44 


IM Vir 


F8 


0.67 ± 0, 


.04 


1949 ± 578 


205 ± 130 


0.367 


1.02 ± 0. 


.06 


2026 ± 603 


311 ± 121 


0.64 


27.9 


13.4 




F9 


0.77 ± 0, 


.02 


3474 ± 680 


1323 ± 277 




1.27 ±0. 


,04 


3460 ± 289 


2041 ± 401 




63.7 


61.2 


V471 Tau 


FIO 


1.44 ± 0, 


.04 


2874 ± 605 


595 ± 204 


0.96 


2.12 ± 0. 


,06 


2136 ± 286 


580 ± 109 


1.28 


42.4 


17.4 




Fll 


1.31 ± 0, 


.03 


1272 ± 301 


864 ± 189 




1.70 ± 0. 


,02 


1847 ± 619 


827 ± 307 




38.3 


17.2 


CC Eri 


F12 


4.9 ± 0, 


.3 


666 ± 162 


256 ± 100 


3.86 


4.6 ± 0. 


,2 


781 ± 98 


336 ± 62 


2.36 


2.0 


0.34 




F13 


6.9 ± 0, 


.2 


3027 ± 441 


493 ± 60 




6.4 ± 0, 


,3 


2844 ± 312 


722 ± 61 




6.7 


2.82 




F14 


7.9 ± 0, 


.1 


4800 ±188 


677 ± 91 




7.4 ± 0. 


,2 


6992 ± 231 


916 ± 61 




6.7 


6.26 


BP Eri 


F15 


1.59 ± 0, 


.04 


939 ± 340 


897 ± 439 


1.20 


2.48 ± 0, 


,10 


1094 ± 293 


967 ± 201 


1.90 


1.4 


0.56 




F16 


1.72 ± 


,04 


4476 ± 867 


2800 ± 1103 




2.54 ± 


()4 


4062 ± 614 


3200 ± 780 




3.1 


2.76 




F17 


1.68 ± 0, 


.03 




820 ± 136 




2.61 ± 0, 


,03 




1014 ± 697 




0.8 





Note: FN is flare name, Aq is count rate at flare peak, is flare decay time, T-r is flare rise time, q is quiescent state count rate, Lx is total 
X-ray luminosity during the flare and Etot is total X-ray energy emitted during the flare 
*Lx for each flare is integrated using the values given in Table 6 based on spectral fits 



00 



Table 5. Parameters obtained from the fitting of flare light curves observed with PN detector in the soft, 
medium and the hard bands. 



Object 



Soft (0.3-0.8 koV) 



Medium (0.8-1.6 keV) 



Hard (1.6-10.0 keV) 



s 



Name 




(-! 




{'■■ 


nil. 


Ao 


L J 


q 




(■-) 


{'■■ 


.UUI.- s 


1 J 


q 


>.-) 


Tj- 

'■>■': 


>.•> 


Ao 




q 


a 


V368 Cep 


Fl 


1486 ± 171 




3 


.96 


± 0. 


.09 


2.69 


1106 ± 110 




3. 


.64 ± 0, 


.09 


2.15 


840 ± 97 




0, 


.74 ± 0. 


.03 


0.31 


o 

e 




F2 


2749 ± 747 


374 ± 79 


3, 


.90 


±0. 


.07 




1949 ± 318 


249 ± 43 


3. 


.86 ± 0. 


.07 




1546 ± 212 


186 ± 37 


1 


.35 ± 0. 


.06 








F3 


318 ± 47 


545 ± 90 


4, 


.21 


± 0. 


.10 




610 ± 113 


707 ± 184 


3 


.46 ± 0, 


.10 




324 ± 46 


283 ± 56 


0, 


.81 ± 0. 


.07 








F4 


1258 ± 264 


545 ± 90 


3, 


.46 


± 0. 


.10 




1184 ± 166 


792 ± 201 


3 


.02 ± 0, 


.07 




1130 ± 330 


709 ± 224 


0, 


.52 ± 0. 


.04 




1 


XI Boo 


F6 


6662 ± 748 


2821 ± 236 


2, 


.62 


± 0. 


.04 


1.40 


4666 ± 340 


3429 ± 291 


3 


.74 ± 0. 


.07 


1.22 


2120 ± 191 




0, 


.47 ± 0. 


.02 


0.06 




F7 


lOlBB ± 786 


4120 ± 643 


2, 


.66 


± 0. 


.03 




10599 ± 543 


3320 ± 296 


3 


.86 ± 0. 


.04 




6772 ± 492 


2683 ± 556 


0, 


.46 ± 0. 


.01 




IM Vir 


F9 


3810 ± 466 


1194 ± 296 


0, 


.56 


± 0. 


.02 


0.32 


3226 ± 427 


848 ± 112 


0. 


.56 ± 0, 


.02 


0.27 


790 ± 130 


360 ± 171 


0, 


.21 ± 0. 


.01 


0.06 




V471 Tau 


FIO 


2946 ± 530 


861 ± 146 


0, 


.98 


± 0. 


.03 


0.66 


2466 ± 464 


676 ± 179 


0, 


.98 ± 0, 


.03 


0.53 


1409 ± 549 


806 ±261 


0, 


.18 ± 0. 


.02 


0.07 






Fll 


1962 ± 267 


1203 ± 272 


0, 


.87 


± 0. 


.03 




1264 ± 268 


1149 ± 209 


0. 


.77 ± 0, 


.03 




278 ± 121 


247 ± 106 


0, 


.16 ± 0. 


.02 






CO Bri 


F12 


591 ± 149 


383 ± 99 


1 


.62 


± 0. 


.08 


1.27 


641 ± 90 


503 ± 151 


2. 


.22 ± 0, 


.08 


1.27 


452 ± 138 


460 ± 123 


0, 


.33 ± 0. 


.04 


0.07 






F13 


3182 ± 501 


800 ± 115 


2 


.99 


± 0. 


.09 




2429 ± 232 


712 ± 92 


4. 


.23 ± 0, 


.05 




2246 ± 661 


692 ± 163 


1 


.03 ± 0. 


.10 








F14 


6446 ± 334 


1662 ± 138 


3, 


.40 


± 0. 


.06 




6765 ± 360 


1262 ± 194 


3. 


.87 ± 0, 


.10 




3934 ± 262 


513 ± 93 


1 


.14 ± 0. 


.06 






EP Eri 


F16 


1223 ± 256 


609 ± 194 


1 


.61 


± 0. 


.03 


1.21 


715 ± 266 


601 ± 201 


0. 


.91 ± 0, 


.07 


0.68 

















Note: FN is flare Name, Aq is count at flare peak, r^^is flare decay time, Tr-is flare rise time, q is the quiescent state count rates. 



A study of X-ray flares - I. Active late-type dwarfs 29 

Table 6. Spectral parameters derived for the quiescent emission of each program star 
from the analysis of PN spectra (MOS spectra for XI Boo) accumulated during the time 
interval 'Q' as shown in Fig. [1], using APEC 2T model. 



Object N^^ Z kTi EMi kTs EM2 Lx xl (DOF) 



V368 Cep 


1 q+0.4 


14+001 
'-'•-'-^-0.01 


S5+°-°^ 

'J-'-'^-O.Ol 


7 4+0.5 


90+°-°^ 


K n+0.2 
'-'•-'-0.3 


6.8 


1.71(429) 


XI Boo 


< 0.95 


< 0.23 


20+0-08 


1 CI +0.70 
-'-•^-'--0.51 


fi7+o.oi 
u-J' -0.01 


1 2+01 


0.51 


1.68(72) 


IM Vir 




'-'•-'^^-0.02 


40+0.05 


7 4+1-2 
' -^-1.1 


0.971H^ 


0+1.0 

8.d_23 


8.3 


1.13(280) 


V471 Tau 






41 

^•"^^--0.02 




0.90^°:°^ 


fi ci+0-4 
0-0-0.4 


7.5 


1.16(341) 


CC Eri 


1 0+0.8 
-'-•'-'-0.7 


n 1 0+0.02 


SI +0-01 

U-'JJ--0.01 




0.83in? 


5+0-2 


0.9 


1.05(302) 


EP Eri 


< 0.4 


"-'•^-'--0.01 


Sfi+0-01 




0.69in^ 


3+0-2 


0.7 


1.44(291) 



Note: Nh is in 10^° cm~^ temperatures (kT) are in keV, emission measures (EM) 
are in 10^^ cm~^, and X-ray luminosity (Lx) is in 10^^ ergs s~^. xl is the minimum 
reduced and DOF stands for degrees of freedom. 
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Table 7. Spectral parameters for each temporal segment of the flare in the targeted stars. 
The temporal segments of the flare are marked by vertical lines in Figure [H 



Object LS FS kT EM Lx xl (DOF) 



V368 Cep Fl Dl 1.201°;^^ 3.71°;^ 9.4 0.97(213) 

D2 IMt^M l-Q^as 8-0 1-11(197) 

D3 0.74t°:23 0.9l[!J 7.2 0.94(186) 

D4 0.58^°:?^ 0.51°:^ 6.7 0.82(170) 

F2 Rl 2.40i[Ji^ 3.0l[!:4 9.7 1.19(157) 

D5 2.6lt[|:| 4.7t[ji 11-8 1-14(220) 

D6 2.46;°;2^ 3.5l°;^ 9.5 1.02(163) 

D7 1.86tS;43 3.41°:^ 8.8 1.12(163) 

F3 .. 0.98t[|:|]^ 4.6l|]J 9.5 1.18(234) 

F4 R2 1.79+°i^ 2.6+°;^ 8.2 1.06(255) 

D8 1.22+°:°? 2.5+J5J 8.2 0.88(183) 

D9 0.7lt[!:J? 1.6l[!J 7.1 1.04(174) 

Dio o.57^[J:}^ i.ol[!:^ 6.4 0.96(220) 

XI Boo F5 .. 0.8i;[|:^^ 0.10^[J:[J^ 0.54 1.37(89) 

F6 Rl 0.94t°:°^ 0.34t°:°2 0.74 1.17(103) 

Dl 0.931°:°^ 0.68+°:°^ 1.02 0.97(72) 

D2 o.8otE5:J54 o.57+E5:|j4 0.93 0.90(65) 

D3 0.45l|]:|J^ 0.83 0.87(61) 

D4 0.78+°:°^ 0.39^°;°^ 0.78 1.21(60) 

F7 *R2 1.08;U;?^ 0.28;[|:J^ 0.93 1.27(93) 

*D5 1.15t°:}^ 0.42+°:^^ 1.11 1.20(85) 

*D6 1.2ll°:i^ 0.43t°:°^ 1.07 0.93(82) 

*D7 1.04+°:^^ 0.33+[J:^^ 0.99 0.95(81) 

*D8 0.95+[J:J]^ 0.34^°:°^ 0.88 1.09(78) 

D9 0.77^[J:[!^ 0.29;EJ;|]^ 0.73 1.07(73) 

IM Vir F8 Rl < 5.16 7.81?;^ 14.6 1.10(60) 
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Table 7 (cont'd) 



Object LS FS kT EM Lx xl (DOF) 



CC Eri 





Dl 


q nq+i oD 
"^■^"-•-0.91 


Li 


13.3 


1.23(88) 


F9 


R2 






15.9 


1.08(103) 




D2 


, , _i_n 1 n 
1 oq-|-U.iU 




15.8 


0.86(36) 




D3 


0.891°:}^ 


p- O+U.S 

'J-O-0.8 


11.8 


1.22(77) 




D4 


l-02lgi3' 


^ ^ _i_n T 

3.01°:? 


10.3 


1.12(81) 




D5 


o-72j:°:i^ 


2.81°:? 


9.90 


0.79(31) 


FIO 


Rl 


-'-■'^"^-0.12 


4.4i°:« 


10.6 


0.85(49) 




Dl 


^ _ _ 1 n 1 /I 

i-29tg:l^ 


5.61°:? 


12.3 


1.18(74) 




D2 


1.24^^4 


q c+o.e 
"J-^-o.e 


10.3 


1.03(64) 




D3 


i.08+°:?2' 


9 q+U.O 
^•"-•-D.S 


9.2 


0.80(80) 


Fll 


R2 


_l_2 


„ in 


10.8 


1.14(43) 




D4 


0-911^:?^ 




9.7 


0.89(50) 




D5 


< 2.79 


1-31^:^ 


9.5 


1.35(56) 




D6 


< 1.48 


o.rl[i:i 


8.3 


0.79(72) 


F12 


Rl 


> 1.03 


n 1 9+<J-Uo 

^•-■^^-0.05 


1.0 


0.85(74) 




Dl 




^ -. . -1-0 Ci^ 

'-'•J-^-0.04 


1.0 


1.24(135) 


F13 


R2 


1 n 1 o 

1 oq+u.i^ 

-'-■^'5-0.12 


n 1x9+0.05 
'-'•^'^-0.05 


1.3 


1.14(129) 




D2 


-1 17+0.11 
^■^'-0.18 


r, 77+0.07 
'-0.07 


1.5 


1.24(82) 




D3 


1 9q+0.25 




1.3 


1.26(86) 




D4 


1 no"^"-^'^ 

■'-•'-''-'-o.io 


^•^^-0.05 


1.3 


0.87(101) 




D5 


^■^ ' -0.07 


4q+o.05 

'-'•^'-'-0.05 


1.2 


0.91(126) 


F14 


R3 


-r 90+0.12 
^-•^'-'-0.12 


0.521°:°^ 


1.3 


1.14(129) 




D6 


1 90+0.35 


0.801°:°^ 


1.5 


1.32(120) 




D7 


-1 9-1+0.07 
-'-■^-'--0.07 


0.701°:°! 


1.4 


1.20(218) 




D8 




0.591°:°! 


1.3 


1.27(209) 




D9 


u.yo_QQg 


0.311°:°! 


1.2 


1.17(217) 
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Table 7 (cont'd) 



Object LS FS kT EM Lx xl (DOF) 



EP Eri F15 


Rl 




n 1 3+°°^ 
'-'■^'J-O.OS 


0.7 


1.24(74 ) 




Dl 




A -1 7+0.09 
-0.05 


0.7 


1.31(126) 


F16 


R2 


0.95i°:}° 


A OO+0.04 
U-28_o.04 


0.8 


1.14(185) 




D2 


0.991°:}J 


A 97+0.02 
U-^ ' -0.02 


0.8 


1.10(187) 




D3 


0.781H.^ 




0.8 


1.16(195) 




D4 


' -^-0.08 




0.7 


1.01(140) 


F17 




rv o -1+0.08 
U-O-l— 0.08 




0.8 


1.05(190) 



Note: LS stands for light curve segments, FS stands 
for "Flare Segments" , Nh is in units of 10^° cm~^, tem- 
peratures (kT) are in keV, emission measures (EM) 
are in 10^^ cm"^, and X-ray luminosity (Lx) is in 10^^ 
ergs s~i. xl is the minimum reduced x^ and DOF 
stands for degrees of freedom. 

*Spectral parameters of hot component in APEC 2T 
fit. 
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Table 8. Parameters derived for flares. 



Object(FN) C T,nax L p n* V Ef, B /3 NL 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 



V368Cep(Fl) 


0.72 


± 


0.03 


2, 


.66 


± 


0.25 


6.0 ± 0.6 


1.02 


± 


0.32 


1.44 


± 


0.20 


1.78 


2.36 


160 






V368Cop(F2) 


0.80 


± 


0.40 


6 


,20 


± 


0.80 


19.8 ± 8.6 


4.37 


± 


2.55 


2.56 


± 


1.16 


0.72 


4.75 


331 


0.12 


~ 2 


V368Cep(F4) 


1.60 


± 


0.20 


4, 


,03 


± 


0.18 


11.7 ± 1.3 


2.05 


± 


0.36 


1.84 


± 


0.32 


0.74 


3.06 


227 


0.27 


^ 8 


XIBoo(F6) 


0.73 


± 


0.30 


1. 


.89 


± 


0.08 


18.6 ± 7.6 


0.13 


± 


0.06 


0.25 


± 


0.10 


10.54 


0.08 


58 


0.61 


~ 26 


XIBoo(F7) 


2.1 


± 


0.4 


2 


53 


± 


0.27 


79.4 ±6.1 


0.07 


± 


0.02 


0.11 


± 


0.02 


37.81 


0.01 


43 


0.11 


~ 1 


IMVir(F9) 


0.60 


± 


0.20 


4, 


,66 


± 


0.37 


18.3 ± 8.9 


2.02 


± 


1.10 


1.57 


± 


0.78 


4.55 


2.07 


226 


0.34 


~ 12 


V471TauF(10) 


1.16 


± 


0.07 


3. 


,36 


± 


0.33 


21.7 ± 2.6 


0.63 


± 


0.20 


0.68 


± 


0.10 


8.88 


0.46 


126 


0.37 


~ 14 


CCEriF(13) 


0.51 


± 


0.25 


2, 


,57 


± 


0.29 


< 14 


> 


0. 


43 


> 


0. 


60 


< 2.15 


> 0.41 


103 


0.33 


~ 11 


CCEriF(14) 


0.,54 


± 


0.15 


2 


70 


± 


0.21 


20.0 ± 10.0 


0.35 


± 


0.20 


0.47 


± 


0.25 


3.60 


0.26 


95 


0.26 


~ 7 


EI'Ei iffii,, 


i.id 




(|.:!IJ 




112 






2:',.r^ = l.r, 






li.or, 


11. L!:^ 




(|.IJ,''> 




(1.(17 






7 



(1) FN is flare name 

(2) C is slope of decay path in density-temperature diagram, 

(3) Tmax is the maximum temperature in the loop at the flare peak in unit of 10^ K based on spectral flt and equation (5), 

(4) L is half-length of the flaring loop in unit of 10® cm and is determined using equation (3), 

(5) p is the maximum pressure in the loop at the flare peak in units of 10^ dyne cm~^ (sec equation (7)), 

(6) Tie is the maximum electron density in the loop at the flare peak in unit of lO"''"'' cm~^, assuming that hydrogen plasma is totally 
ionised (p = 2nekTmax)j 

(7) V is volume of flaring plasma in unit of lO'*^ cm"* and estimated using the equation (6), 

(8) Ef^ is heating rate per unit volume at the flare peak in the units ergs s""'" cm"'^ estimate from the loop scaling laws {see equation 
(8)). 

(9) B is minimum magnetic field necessary for confinement in Gauss. 

(10) /3 is loop aspect ratio (r/L). 

(11) NL is number of loops involve in the flaring plasma. 

(*) The pressure as derived from equation (7) is a maximum value, appropriate for the equilibrium condition, which can be reached 
only for very long-lasting heating. The density and volume values derived thereby are therefore be treated as upper and lower limits, 
respectively. 
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Figure 1. MOS and PN Light curves of six G-K dwarfs 
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Time (in s since JD 24532 18.78985) 
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Fig. [T] Continued 
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V368 Cep (M0S1+M0S2; Bin lime 200.0 s) 
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Figure 2. (a) MOS and PN light curves at three bands soft (0.3-0.8 keV), medium (0.8-1.6 keV), and hard(l. 6-10.0 kcV) and 
hardness ratio HRl and HR,2 curve, where HR,l=medium/soft and HR2=hard/medium of V368 Cep. 
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Fig.EIb) - Similar to Fig. [2] (a), but for XI Boo (only MOS) 
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IM Vir (M0S1 +M0S2; Bin time 200.0 s) IM Vir (PN; Bin lime 200,0 s) 
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Fig. [He) - Similar to Fig. [2] (a), but for IM Vir 



V471 Tau (M0S1+M0S2, Bin time 200s) V471 Tau (PN; Bin Time 200.0 s) 
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Fig. Eld) - Similar to Fig. [2 (a), but for V471 Tau 
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Figure 3. Colour-colour diagram (HRl versus HR2) for the various time segments (Q: quiescent state, Ri{i=l,2..): rise phase, 
Di(i=l,2..): decay phase) observed for each source. The curves overplotted on the data are for model simulation where two 
temperature plasma is assumed to predict the colors (see i|4.2l for details). 
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Figure 4. Quiescent state spectra of PN data with APEC 2T model in upper subpanels of each graph, contributions are 
given in the lower subpanel of the each graph. 



42 J. C. Pandey and K. P. Singh 




Figure 6. Plots of (a) temperature and (b) emission measure of dominant component of 2T fit during the flares F6 and F7 of 
the star XI Boo. 



A study of X-ray flares - I. Active late-type dwarfs 43 





Figure 8. Plots of (a) temperature and (b) emission measure during the flares FIO and Fll of the star V471 Tau. Here, arrow 
shows the upper limit on data point. 
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Figure 11. The density-temperature diagram, where EM^/'^ has been used as a proxy of density. Symbols for the each 
flares are given at the bottom right corner. Straight lines represent the best linear fit to the corresponding data, f is slope of 
density-temperature diagram. 
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Figure 12. Plot of peak flare luminosity (Lx/) versus stellar quiescent luminosity {Lxq) 



